On Spherically Symmetric Motions of a Gaseous 
Star Governed by the Euler-Poisson Equations 

Tetu Makino 
December 31, 2012 



1 Introduction 
1.1 

We consider spherically symmetric motions of a gaseous star governed by the 
Euler-Poisson equations: 



dp dp du 2 

wr + u— + p— + -pu = 0, 

at or or r 

,du du. dP 9$ ,„ , 
p{ m +U Tr ) + ^ = - p ^ <?<*><> <r)> 

Here p is the density, u the velocity, P the pressure, $ the gravitational potential, 
and go is the gravitational constant. In this work we assume 

P = Ap\ (2) 

where A and 7 are positive constants, and we assume 1 < 7 < 2. 

Introducing the mass 

m := 47r / p(t,r')r' 2 dr' , 
^0 

we can write the equations as 



On the other hand, equilibria for the equations (1) are governed by the 
ordinary differential equation 

-L—(——\ - 4 ng p 

r 2 dr \ p dr J 

which is called the Lane-Emden equation. The following result on this Lane- 
Emden equation is well-known: 

Lemma 1 Assume 6/5 < 7 < 2. For any positive number po given, there is an 
equilibrium p — p(r) with positive numbers R, p\ such that p{r) is positive and 
analytic in < r < R and 

p(r) = p (l + [r 2 ]i) asr^O, 

p(r) = Pl {R-r)^{\ + [R-r, (R-r)^}^ asr^R-0. 

Notational Remark. Here and hereafter [X] q denotes a power series of the 
form J2j> q a jX-' with positive radius of convergence, and [X,Y] 9 a convergent 
power series of the form J2j+k>q a jkX^Y k . 

Remark. In the expansion of p(r) as r — > R, the terms including 

(R — r) i- 1 actually appear if is not an integer. Let us prove it. Otherwise 

7-1 

we would have 

p(r) = pi( J R-r)^r(l + [ J R-r]i) 

and the function 

U(r) := pW 1 = P r\R -r)(l+[R- r] x ) 

would be analytic at r = R. Now U satisfies 

d 2 U 2 dU _ 47r. g o(7 - 1) 
dr 2 r dr 

Since U is analytic, the left-hand side is analytic, and so, the right-hand side 



A-y 



-p 1 (R-r)y=*(l+[R-r] 1 ) 



would be analytic at r = R. Then shoud be an integer. This contradicts 

7- 1 

7 1 

to that = h 1 is not an integer. 

7 — 1 7 — 1 

Moreover if we write R = R(po),pi = pi(po) as functions of po, then we 
have 



i?( Po ) = Po 2 RQ), pM = p^- 1} P i(i). 

2 

In fact p(r) = K-<- 2 p(r/K) is an equilibrium, too, for any K > and any 
equilibrium p(-). 

For a proof of Lemma 1, see, e.g., [4], and [5], Chapter V or [9 , Chapter IX. 
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1.2 

In the following discussion we assume that 6/5 < 7 < 2 and we fix an equilib- 
rium p(r) with the properties in the above Lemma. Put P(r) = Ap(r) 7 . 

We are going to construct solutions around this fixed equilibrium. To do so, 
we take the Lagrangian mass coordinate m as the independent variable instead 
of r. Then we can write the equations as 

^+4wp 2 {r 2 u) m =0, 



9U . On 711 

^+4ttt P m = -g - 

,47T,/ P ' 



V47T ./n > 

Since 

dr dr 1 

u, 



dt dm Airpr 2 ' 

the equations are reduced to the single second order equation 



r tt +4irr 2 P m = -g ^, (4) 



where 

P = AUnr 2 pY\ 
\ dm J 



Now we derive the equation for the perturbation y defined by 

r = f(l + y). (5) 

Here f (to) is the function of the Lagrangian mass variable m associated with 
the fixed equilibrium. In other words, it is the inverse function of 



f P{r') 
Jo 



Keeping in mind 



we have 



m = m(r) = An p(r )r dr' . 



dr dr / r dy \ 

dm dm\ ^ f m dm)' 



Here G(y, v) = 3-yy + jv + [y, v] 2 is defined by 

(l + y)-^(l + y + v)-i = 1-G(y,v). 
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Then the equation is reduced to 
where we have used 



d did 



dm dr Airpr 2 dr 
We note that the equilibrium satisfies 

IdP m 
p dr r z 

Let us introduce H(y) = Ay + [y] 2 by 



H(y) = (l + y) 



2 



{l+yf 

Then the equation can be written as 

Here we have ued the abreviations r, p, P, P r instead of f, p, P, dP/df. We 
consider this nonlinear wave equation. 

It is easy to verify by a scale transformation of variables that we can assume 

R = l, A=- 

7 

without loss of generality. Hence we shall assume so, and we consider the 
equation (6) on the interval (0, 1) of the variable r, and we assume that 

p = p (l + [?" 2 ]i) asr^O, 

1 7 

p = pi(l — r)~(l + [1 — r, (1 — rj^ji) as r -> 1 

and 

P= ip 7 . 

7 

2 Analysis of the linearized equation 
2.1 

The linearized equation is 
d 2 y 

d£+£y = °> ( 7 ) 

4/ := -^W 3 ^ + ^|)) + ^-^ 



pr dr \ \ ' dr/ ) pr 
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and the associated eigenvalue problem is 

Cy = Xy. 

Let us use the Liouville transformation: 



^ ~ Jo V-Tp^ 7 *' V-=r 2 (jPp)^y. 
Through this transformation the equation 

Cy = Xy + f 



turns out to be the standard form 
where 

7 P/2 / 7-3 7 l+jrm r \ Pr (7 + 1)(3 - 7) ( Pr\ 2 ^ 

q p \r 2 V 2 4 m ) rp 16 V p ) J 

g = r 2 (jPp)if. 



The variable £ runs on the interval (0,£ + ), where 
Since 



i+ ■= [ \/-^dr < 00. 



— — r as r — > 0, 

7-Po 



we see 

Since 
P 

and 

we see 



7P0 2 2 
<7 o "72 as £ ->• 0. 



— (1-r)- 1 , ^~p7 _1 (l-r) asr^l-0, 

- 1 P 



l-r-^" 1 ^-^) 2 as£^£ + -0, 



7 P(7 + l)(3-7)/Pr\ 2 l(l+7)(3-7) 1 



* p i6 v P y 4 (7-1)2 ( e + - c) 2 

as £ -> £ + - 0. It follows from 1 < 7 < 2 that 

1 (1+7X3-7) 3 

4 (7-I) 2 4' 
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Of course q is bounded from below, but it is difficult to know whether its 
minimum is positive or not. Anyway, the both boundary points £ = 0, are 
of limit point type. See, e.g., [5], p. 159, Theorem X.10. The exceptional case 
7 = 2 will be dicussed later. Hence we have the following conclusion: 

Proposition 1 The operator T ,T>(T ) = Cg°(0, T r] = —r]^ + qr], in 
L 2 (0,i^+) has the Friedrichs extension T, a self-adjoint operator, whose spec- 
trum consists of simple eigenvalues \\ < ■ ■ ■ < X n < X n +i <•••—> +oo. 

In other words, the operator Aq, T>(Aq) = C§°(0, 1), A Q y = Cy in L 2 ((0, l),r 4 pdr) 
has the Friedrichs extension A, a self-adjoint operator with eigenvalues (X n ) n . 

The domain T>(T) of the Firedrichs extension T is, by definition, 

P(T) = {?? e L 2 ((U+) I 30„eC o °°(O,e + ), Q[0m-0»]^O 

as Tn, n — ^ oo, <p n — > r\ inL 2 (0,£_|_) 
and — rj^ + qr) G L 2 (0, £+) in distribution sense}, 

where 

Q[4>] := / + f # 2 + (g + c)|0| 2 Ve, 



and c is a constant > | min<7|. But T>(T) is characterized as follows: 

2?(T) = {77 e C[0,£+] | r?(0) = V (t+) = 0, -r fe + 977 G L 2 (0, £+)}. 

Let us prove it, denoting by M the right-hand side. Let 77 £ T>(T). Then 
there are </>„ G Co°(0,£+) such that 4> n — > r\ in L 2 and Q[0 m — 0„] — > 0. Since 

|^n(0 " K{t)\ < \Tt I {{<t>m - <Me?d£ < JlQ[<t> m - <t>n] -> 0, 



we have <pn V uniformly on [0, £+]. Hence 77 € C[0, £+] and 77(0) = 0. Similarly 
r)(( + ) = 0. Thus V{T) C M. Let r/ G A/. Put / := -77^ + 977 G L 2 . Then 
—77^ + (g + c)?7 = g := f + en G L 2 . since belongs to the resolvent set of 
T + c, we have u := (T + c) _1 s G P(T). Hence w := 77 - v G C[0,£+] and 
jyj(0) = ttj(£ + ) = 0, -io^ + (q + c)w = 0, for X>(T) C M. Using q + c > 0, we 
can deduce that w = and 77 = u G T>(T), that is, M C V(T). (In fact, if w did 
not vanish identically, there would exist a G (0,£+) such that Dw(a) = and 
w(a) 7^ 0. If 7ij(a) > 0, then 

Dw(0= [ D 2 w(e)d?= [ (q + c)w(e)dt' 

J a J a 

implies Dw(£,) > for a < £ < £+ and it contradicts to 7-i>(£+) = 0. If 77j(a) < 0, 
then 
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implies Dw(£) > for < £ < a and it contradicts to w(0) = 0.) ■ 

Although it is not easy to judge the signature of min q, we have the stability 
result by S.-S. Lin, [B], 1997, as follows: 

Proposition 2 If and only if 4/3 < 7 < 2, the least eigenvalue Ai is positive. 
Proof: The function y = 1 satisfies 

Cy = — (4-3 7 )P r >0. 
pr 

Let us consider the corresponding function 

Tli = r {jPp) 4 

thorugh the Liouville transformation. It is easy to show that rji and dr/i/d^ 
vanish at £ = 0, £+ and 771 € D(T). Let 0i(£) be the eigenfunction of — d 2 /d^ 2 + 
Q associated with the least eigenvalue Ai. We can assume </>i(£) > for < 
£ < £+ and </>i and dcj>i/d£ vanish at £ = 0,£+. Then the integration by parts 
gives 

ri+ /■£+ 
Ai / 0i»7id£ = / + OTi)*;- 

Jo Jo 

Since 

-ViM + 91i = r(7F)3p~4(4 _ 3 7 )p r 
and P r < 0, we have the assertion. ■ 



2.2 



Let us introduce the variable x defined by 
tan 2 9 



1 + tan 2 9 



2 



sdr, 



where k — ?r/£+. Then x runs over the interval [0, 1] while r runs over [0, 1], 
and 

dx 1— r -T+i 



Since 



dr 



— Ky/x(l — x)p 



(J_ 

dr 
d 2 
dr 2 



Kyjx{l - xjp 2 



-,+1 d 2 



— = K X(l - x)p 2 —- 

+ (± K 2 (l-2x)p-^ 1 



-7 + 1 



ci/a;(l - x)p 2 p r \ — 



dx 
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we have 

k Cy = -x(l-x)— + 

/I/ n 7+11 i— 7-3 4 1 i— ~,-i\dy 



As r — > 0(a; — > 0) we have 



Then it follows that 



2 7- 1 

-p 2 V^ 1 + Mi)> 



'^■^/xlT^xjp 2 ^ = 2V1 - s(l + [a;]i). 



K - 2 



£y = - a! (l + 1^)1 + Hoy. 



On the other hand, as r — > l(x — > 1), we have 

1-x = KV 7+1 (l-r)(l + [l-r,(l-r)^r]i), 

1-r = K- 2 pJ- 1 (l-x)(l + [l-a:,(l-a;)^T]i), 

= — ^-(l-r)^(l + [l-r,(l-r)^H). 
7-1 

Then it follows that 

n~ 2 Cy = —x(l — x)^ r ^ + 
ax"' 

+ (--__y + I 1 - C 1 - aO^h) ^ + I 1 - x, (1 - ar)^] ». 

Changing the scale of t, we can and shall assume that k = 1 without loss of 
generality. 

Summing up, we have: 
Proposition 3 We can write 

Cy = -x(l - ,)g + (§(1 -x)-£.)| + + L (^, 
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where 

Lx(x) = 

N 

L (x) = [x] , 

= [1 - a;, (1 -a;) ^ Jo- 
Here N is the parameter defined by 

N=^j * 7 = 1 + ^. (9) 

2.3 

Let us consider the simplified linear operator 

d 2 /5„ , iV \ d 



A := a; 



P-^ + ^-^i^ (10) 

= x 2 (l-x) 2 + — a; 2 (1 - x) 2 — 
dx dx 

The function y = 1 is a solution of Ay = and another independent solution of 
Ay = is given by 



y = ij)(x):= I X~i(l- Xy^dX. 



We take the Hilbcrt space X of functions of < x < 1 defined by 
X = {y | \\y\\ x <oo}, 

1/2 



IHIx = ( J IvixtfxHl-x^dx)' 

Then y = 1 belongs to X but ^ does not belong to X. It is easy to prove directly 
the following 

Proposition 4 Let y e X satisfy Ay = f g X. T/ien t/ie following formulae 
hold: 

y(x) = i,{x) I* f(X)xi(l-X)%- 1 dX + 
Jo 

+ [ ^(X)f(X)xl(l-X)^- 1 dX + C, (11) 

J a: 

/ ./(X)X2(l-X)^- 1 dX = 0, (12) 
Jo 

a;-*(l-a:)-T [* f(X)Xi(l- xfi^dX (13) 
Jo 

-a;-5(l-X)-^ A /(X)xi(l -X)-^" 1 ^. (14) 

J X 



dx 



Here C is an arbitrary constant. 
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2.4 

Now let us fix a positive eigenvalue A = A„ and an associated eigenfunction $(r) 
of£. Then 

2/i(*,r) = sin\/A(i-*o)$(r) 
is a time-periodic solution of the linearized problem. 
Moreover we can claim that 

$(r) = C Q (l + [r 2 ]!) asr^O, 

and 

$(r) = Ci(l + [1 - r, (1 - r)^r]i) as r 1 - 0. 

Here Co and C\ are non-zero constants. Other independent solutions of Cy = Xy 
do not belong to L 2 {r A pdr) both at r ~ 0,r ~ 1. To prove this, we use the 
following lemma: 

Lemma 2 Let ms consider the equation 

6(2, z a ) = a + [z, z°]i, c(z, z a ) = [z, z a ] , 

and let the positive number a satisfy a > 2. Then 1) there is a solution yi of 
the form 

y x = l + [z,z a ]i, 
and 2) there is a solution yi of the form 

y 2 = z- a+1 (l + [z,z a ] 1 ) 

provided a £ N, or 

y 2 = z- a+1 (l + [z,z a ] 1 ) + % 1 logz 
provided a £ N. Here h is a constant which can vanish in some cases. 
For a proof, see [T], Chapter 4. 

Let us consider the behavior of solutions at 

z := 1 — r -> 0. 

The the equation reads 

A , / 4 , P r .d|/ , (4-37) Pr , P 



f 4 ^ -ay 14-jTj n p 



dz 2 U-z P ; dz (l-z) 7 P a 7P 1 
and 

§ = -ttztzO- + z ^)^ ? = pF 1 ^ + [*, * B ]i)> 

r 7 1 z p 

with a = . Since 1 < 7 < 2, we have a > 2. Hence an application of the 

7-1 

above lemma gives the assertion. Even if N — 4(7 = 2), ?/2 ~ z~^~ +1 does not 
belong to X, and the boundary point r = 1 is of the limit point type. 
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3 Statement of the main result 



We rewrite the equation (6) by using the linearized operator C defined by (8) 

as 



H+( 1 + ^(y,rg))£y + G„(r,y,r^)=0 > 



dy 



dy" 



where 



Gjfav) = (l + y) 2 (l + ^d v G 2 (y,v)) - 1, 

P P r 
Gn(r, y, v) = — ~G II0 (y, v) + —G H1 (y, v), 
pr z pr 

G II0 (y,v) = (l + y) 2 (3d v G 2 -d y G 2 )v 

= -2 7 (l + y)-^+ 1 (l + y + U )-^ 1 « 2 , 

G in (y,v) = + 8 V G 2 • ((-4 + 3 7 )y + 1V ) + 
7 

+H-4y(l + y) 2 -(l + y) 2 G 2 . 



(15) 



Here 



G 2 {y,v) 



G(y, v) - (3jy + jv) = [y, v} 2 , 



d v G 2 := —G 2 = — - 1 =ly,v} 1 . 

We have fixed a solution yi of the linearized equation y tt + Cy = 0, and we 
seek a solution y of (6) or (15) of the form 

y — eyi + ew, 

where e is a small positive parameter. Then the equation which w should satisfy 
turns out to be 

d 2 w / / dw \\ „ , / dw 



ni 2 + (l + ea(t,r,w,r^,ey)Cw + eb(t,r,w,r-^,e) 
= ec(t,r,e), 



(16) 



where 



e 1 Gi(eyi + ew, ev\ + efl), 
(-F/ + f jj) 



a(f, r, u;, 0,e) = 
r, w, fi, e) = 

c(i,r,£) = 

Here v\ stands for rdyi jdr and 

Fj := — e~ Gj{ey\ + ew, ev\ + eSY)£yi, 
F n := -e~ 2 Gn(r, ey 1 4- ew,evx + eQ). 



w =n=o 
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Then the main result of this study can be stated as follows: 

T 

Theorem 1 Assume that 6/5 < 7 < 2 and that is an integer. That is, 

we assume that 7 is either 2,3/2,4/3, or 5/4. Then for any given T > there 
is a sufficiently small positive Eq = Sq(T) such that, for < e < Eq, there is a 
solution w of (16) such that w € C°°([0,T] x [0, 1]) and 



sup 

j+k<n 



dl) KcTr) W 



< C n e, 

L~([0,T]x[0,l]) 



sdt/ \dr; 

or a solution y G C°°([0,T] x [0, 1]) of (6) or (15) of the form 

y{t 1 r)^ey 1 {t 1 r)+0{e 2 ), 

or a motion which can be expressed by the Lagrangian coordinates as 

r(t, m) = f(m)(l + eyi(t, f(m)) + 0(e 2 )) 

forO<t<T,0<m<M. 

Remark. The corresponding density distribution p — p(t,r), where r de- 
notes the original Euler coordinate, satisfies 

p(t,r) > for < r < R F (t), p(t,r) = for R F (t) < r, 

where 

R F (t) := r(t,M) = 1 + e sin V\(t - t )$(l) + 0(e 2 ). 
Since y(t, r) is smooth on < r < 1, we have 

p(t,r) = C(t)(R F (t) - r)^r(l + 0(R F (t) - r)) 
asr^ R F (t) — 0. Here C(t) is positive and smooth in t. 

Our task is to find the inverse image s Jp~ 1 (ec) of the nonlinear mapping <p 
defined by 

V(w) ■= -5^ + (1 + ea)Cw + eb. (17) 
ot z 

Note CP(0) = 0. It requires a property of the Frechet derivative of *$: 

D^(w)h = h u + (1 + eai)£h + ea 2 oh + ea 2 irh r , (18) 



where 



ai{t,r) = a\t,r,w,r—,e), 
, . da db 

O20 = ^ ! 

a2l(< ' r) = ^ £w+ 9Ti- 
Here Q is the dummy of rdw/dr. We shall use the following observation: 
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Proposition 5 We have 
where 

dy dY 

Y = y 1 +w, y = eY, v = r—=er—. 

or or 



Proof. Since 

:>2 



da _ dd _ (l + y) 2 a ., r , 

Oil ov 7 

^ = dGl £ | c - ldGn 

we have 

a-21 = d v G 2 C{yi + w) + e — — . 

7 ov 

On the other hand 



dGiio 
dv 

and 



= -2 7 (1 + y)" 27+1 (l + y + v)-~<- 2 • (2(1 + y) + (-7 + l)v)v 
dG in _ (1 + V)\ 2 , 



dv 



-diG 2 -{{-A + ^)y + 1 v). 



Thus 



02i = (1 + ^ d*G 2 ■ (CY + —((-4 + 3 7 )y + jV) 
7 V pr 

-^^(1 + y)-^ +1 (l +y + v)-~<- 2 {2(l + y) + (-7 + l)v)v 



pr 

>P„o„ (d 2 Y AdY 



< 1 + y) 2 ^ G2 iwr 



+ 



dr 2 r dr I 

-^f (1 + y)" 27+1 (i + y + v)->- 2 (2(i + y) + (-7 + l)v)V. 

Here V stands for rdY/dr. Using 

d 2 G 2 = -7(7 + i)(i + ?/r 27 (i + y + v)-~<- 2 , 

we get the result. ■ 



4 Proof of the main result 
4.1 

The proof of the main result presented here depends strongly on the preceding 
discussions in [3] and [7]. But a completely parallel treatment might draw dif- 
ficulties even if the linearized operator C is simplified by A, for the structure of 
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A m is complicated. Hence we introduce a cut off function in order to divide the 
difficulties. 

Let us fix a function uj e C°°(R) such that < lj < 1,uj(x) — 1 for x < 1/3, 
and oj(x) = for 2/3 < x. For any function y of < x < 1, we put 

y°(x) =u(x)y(x), y\x) = (1 - u(x))y(x). (19) 

Then, y — y° + y 1 , supp[y°](:= the support of y°) C [0,2/3] and suppfy 1 ] C 
[1/3,1]. 

We consider the operators 

d 2 5 d 

A ° := X d^ + 2Tx (20) 
d 2 N d 

Al ^ Z d^ + YT Z (21) 

Here and hereafter the letter z is used for the variable 

z:=l-x (22) 



First of all we introduce the norm 

\\y\\ m := sup {IK-AoJVIU-.IK-Ai)^ 1 ^-.}. 

0<j'<m 

Clearly we have the equivalence 

Hvllo < Ill/Ik- <2||y||o. 

4.2 

We shall apply the Nash-Moser theorem ([2], p. 171, III. 1.1.1.) to the functional 
spaces 



£ := C°°([0,T] x [0,1]) 



dw 



(B := {we€ | w| t=0 = — =0}, 

at t=o 



and the nonlinear mapping 



= -=-o- + (1 + ea)^w + e6 



Si 

On 2; we introduce the gradings of norms 



(17) 



\\y 



(oo) .. 



sup 

0<a+/3<n,0<t<T 



8 2 \" 



( dt 2 ) 
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and 



where 



0<a+/3<n ' 



3^ \« 

? 



rff 



X,/3 



\\vh,p:={ E (ll(-A )Vlll + IK-aoVIII, 

o<j<p 

Here X = L 2 ([0, l] x , xidx) and Xi = L 2 {[0, l] x , z^^dz). 

Then the Sobolev imbedding theorem shown in 3 implies the equivalence 
of these gradings, say, 

^Nll 2) < < C\\y\C 

where a/2 > max(5/2, N/2). 

Moreover (£ and its closed subspace <£q are tame spaces. In fact, we can 
consider the tame mappings $ : l£ — > (£° x (S 1 : y M> y 1 ) and * : (yo,yi) | -> 
J/o + 2/1) where (£° and (S 1 are copies of the space "<£" in [5] endowed with 
the operators A and Ai respectively. These spaces are tame as claimed in 
[3]. Then €° x (S 1 is tame as the cartesian product of tame spaces ([2], p. 136, 
II. 1.3. 4.) and € is tame as a tame direct summand of a tame space ([2], p. 136. 
II.1.3.3.). 

In order to claim that the mapping *}3 is tame, it is sufficient to use the 
following propositions, which are parallel to |3 ], Propositions 7,8,9: 

Proposition 6 For any meff there is a constant C such that 

dw 



dr 



< C\\w\ 



m+l ■ 



Proof: We have 



dw 
dr 



ry x(l — x)p~ 



. dw 

dx 



2p^x{l + [x\)^ 

i=± , , n dw 
- Pl 2 (l + lz^-l) — . 

Hence [3], Proposition 7 implies the assertion. ■ 

Remark 2: Using the integral representations in Proposition 4, we can 
deduce 

<C||Ay|| i0 o. 

dx 



But it may be difficult to write down a simple integral representation for A m — . 

dx 
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Proposition 7 For any m € N there is a constant C such that 

\\f-g\\ m <C(\\f\Ug\\ + \\f\\ Q \\g\\ m ). 

Proposition 8 If F(x, y) is a smooth function of x and y, then for any m e N 
there is a constant C such that 

\\F(x,y(x)\\ m <C{l + \\y\\ m ). 



4.3 

Let us consider energy estimates for the higher order derivatives of the solution 
of the equation 

D¥(w)h = g, 

where h € <£ , g <E <£, while w € €q is fixed, and 



d 2 h dh 
Dty(w)h = -^jj + (1 + sai)Ch + ea 2 oh + ea 21 r—. 



Using the relations 



uiCh — C{u>h) 
dh d 



, „ dh 
u 2 h + 2wi — , 
or 

rcuih, 



(18) 



where 



4 du 



dr 



1 d 

pr 4 dr 



we can rewrite the equation (18) to the following system of simultaneous equa- 
tions: 

W° = g°-a 1 ^-a h 1 
or 
dh 

Wlih 1 = g 1 +a h°, 
dr 

where 

a 2 d 

9^0 = ttpt + (1 + ea\)C + ea 20 + ea 21 r— + 
dt dr 

d 
or 

d 2 d 
9Ki = ^ + (1 + eai)£ + ea 20 + ea 2 ir— + 

""I Q "0, 

or 

ai = 2(l+eai)cji, 

a = (1 + eai)ui 2 - ea 2 iru> 1 . 
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Moreover we note that 



dx 

N d 

a . \ (i - i,Ao+ r? 



Therefore we can write 



where 



d 2 d 
Tlo = -g-2- o°(t, x)A + 6?(t, + &g(t, x), 



a = (l+eoi)(l-ar), 

,n N ,„ . Li . 1 cte 

61 = _ y + (1 + eai) T + (efl2ir + Q!l) x* : ' 



and 



where 



6q = L + ea 2 o + ao, 



a 2 a 
OTl = ^2 - fll (*> z ) A i + 6 i(*. z ) z ^ + 5 o(*, 2) 



a 1 = (1 +eai)(l - z), 

,1 5 . -Li . . 1 c?z 

&i = ---(1 + eai) V {ea 2 \r - on)- — , 

A Z Z (XT 

bo = L Q + ea 2 o-a . 

First, if |eai| < 1/2, then a > 1/6 for x < 2/3. Since the operator 9Jl 
acts on functions whose supports are contained in < x < 2/3, we can assume 
a° > 1/12 everywhere. 

Since L\ — [x]i, we have Li/x G C°°[0, 1]. 

Since 

r dx 

xTr =2+[x] - 
and since a x = for x < 1/3, we see 6? G C°°([0,T] x [0, 1]). 
Clearly 6° G C°° . 

On the other hand, if |eai| < 1/2, then a 1 > 1/6 for z < 2/3. Hence we can 
assume a 1 > 1/12 everywhere for DJli acting on functions whose supports are in 
< z < 2/3. 

Since L\ = [z, z%]\ and since N/2 e N, we sec Lx/z G C°°([0, l] z ). Although 
r dz =xtl 1 r 1. . 

we know 021 G z • C°° thanks to Proposition 5. By this observation and the fact 
that ai = for z < 1/3, we can claim that b\ G C°°([0,T] x [0, 1]). 



17 



Clearly b\ € C°° . 

Since ao = a± = for x < 1/3 and 2/3 < x (or z < 1/3), there are constants 
A*,B* such that 



dh 1 
1 dr 




< 


A* 


^d/i 1 

Oz 


3Ei 


\\aoh 1 


Wo 


< 


5* 






dh° 
1 <9r 


Xi 


< 


A* 


Ox 




\\a h° 


1*! 


< 


5* 


\h°\\ Xo . 





Then, putting 

A) = ||o°t||i- + VT2|| V^a ^ + &?l|i=° +\/l2A*, 

So - U&8IIL-+B*, 

A x = ||a* t || L =. +^|| N /ia^ + &l|| i oo +VUA*, 

B 1 = \\bl\\ L ~+B*, 

A = m&x{Ao,Ai}, 

B = max{i?o, -Bi}, 

Eo(t) = f ((/i° t ) 2 + a°{t,x)x{h%) 2 )xidx, 
Jo 

Et(t) = f ((h\) 2 +a 1 (t,z)z(h 1 z ) 2 )z%- 1 dz, 
Jo 

E{t) = max{^(t),Si(t)}, 
we have the estimate 

\^ < AE + BE 1 '* J* E(t'y /2 dt> + E^\\g\\ x , . 

Here we recall 

Il9b,o = (||ff°||| + ||9 1 ||| 1 ) 1/2 . 
This is the energy estimate for (h°, h 1 ). 

Similarly the discussion of energy estimates of the higher derivatives of h 
developed in [3] and [7] can be applied to the simultaneous equations for (h°, h 1 ) 
to get the tame estimate of Vty{w) : g n- h. This completes the proof. 

5 Concluding remark 

In order that the equilibrium satisfy that p 7_1 is analytic at the free boundary 
r = R, we assumed that N is even. But 7 = b/3(N = 5) for monoatomic gas, 
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and 7 = 7/5(N — 7) for the air. Therefore it is desired that the result will be 
extended to the case that N is an odd integer. For the case that N is not an 
integer, we may try quite other approach than the present application of the 
Nash-Moser therein. This is an interesting open problem in view of physical 
applications. 
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